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Abstract 

Using the expansions of the heavy meson decay widths in the heavy quark mass 
and QCD sum rules for estimates of corresponding matrix elements, we calculate the 
]j±,o,s an( j g±,o,s meson lifetimes. The results for D mesons are in a reasonable agree- 
ment with the data, while it is predicted: [T(Bd) — T(i? ± )]/rB — 4% (and the lifetime 
difference of the Bd and B s mesons is even smaller); \T(B s s hort ) — T(B l ° n9 )]/T(B s ) ~ 
8% . The role of the weak annihilation and Pauli interference contributions to the life- 
time differences is described in detail. In the course of self-consistent calculations the 
values of many parameters crucial for calculations with charmed and beauty mesons 
are found. In particular, the quark pole masses are: M c ~ 1.65 GeV, M b ~ 5.04 GeV , 
and the decay constants are: f D (M c ) ~ 165 MeV , f B {M b ) ~ 120 MeV . It is 
also shown that the nonfactorizable corrections to the B — B mixing are large, 
Bb{Mi) ~ (1 — 18%) . The values of the unitarity triangle parameters are found 
which are consistent with these results and the data available (except for the NA31 
result for the e'/e which is too large): \V cb \ ~ 4.2 ■ 10~ 2 , \V td \ ~ 1.3 • 10~ 2 , \V ub /V cb \ ~ 
0.10 , { A ~ 0.86 , p ~ -0.40 , r] ~ 0.20 }. 
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1. Introduction 



Experimentally measured properties of heavy mesons offer the possibility to find 
out values of parameters which are of fundamental importance for the Standard 
Model: M c , Mb, V c b, Kt&> Vtd (also fo, Jb-, etc). Below the attempt is described [|TJ of 
self-consistent calculation of all the above parameters using the available experimen- 
tal data on the D and B mesons. The theoretical methods used are: expansions in 
the heavy quark mass for obtaining the effective Lagrangians and QCD sum rules for 
estimates of corresponding matrix elements. The used scheme of calculations looks 
as follows: 

Experiment T(D — > ev + X) : — > determination of M c ; 
M c + Heavy meson mass formulae: — > calculation of M^, 

M c , M b + QCD sum rules: — > calculation of the decay constants fa, fs and the ma- 
trix elements < M\ L e ff \M >; 

Calculation of the D ±,0,s meson lifetimes and predictions for the B ±,0,s meson lifetime 
differences and branching ratios; 

Calculation of the matrix element: < B°\ L™j m9 \B° >, prediction of 
Y(Bf ort ) - Y(B l s on 3). 

Experiment T(B — > ev + X) : — > determination of \V c b\; 
Experiment Am^ : — > determination of \Vtd\] 

Experiment on e K : — > determination of p, rj, V uh , estimates of e'/e, predictions for 
CP-violation in B-decays, etc. 



2. The history 

It is a long standing challenge for theory to calculate the D and B meson decay 
widths. On the qualitative side, two mechanisms were invoked to explain the pattern 
of the D meson lifetime differences: weak annihilation (WA) [[|, and Pauli interference 
(PI) 0,0. As for WA, it was expected that because an admixture of the wave 
function component with an additional gluon or the emission of a perturbative gluon, 
both remove a suppression due to helicity conservation (which leads to Br (n — > 
ev)j 'Br (it — > imp) ~ 10~ 4 ), the D° meson decay width is enhanced. On the other 
hand, it was expected that the destructive PI of two d-quarks (the spectator and 
those from a final state) will suppress the D meson decay width. As for WA, there 
were no reliable calculations at all. For PI, simple minded estimates (see sect. 7) give 
too large an effect which results in a negative decay width. 

For B mesons, it was clear qualitatively that all the above effects which are of 
a pre-asymptotic nature and die off at Mq — > oo, will be less important. However, 
because the patern of the D meson lifetime differences was not really explained and 
well understood, this prevented to obtain reliable estimates of the B meson lifetime 
differences, and only order of magnitude estimates were really available: [5T(B)/T B ] : 
[5T(D)/T D ] ~ 0{f B Mt/p D Ml) ~ O(10^). 
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Moreover, as for WA contributions through perturbative gluon emission (which is 
formally a leading correction ~ 0(Aq C d/Mq) to the deacay width and was expected 
before to be potentially the most important), it has been emphasized recently |J 
that such contributions are of no help at all because, being large (at least formally 
at Mq — > oo) term by term, they cancel completely in the inclusive widths, both 
0(1/ Mq) and 0(1/Mq) terms. It will be shown below (see sects. 8-10) that, never- 
theless, there are important WA contributions but on the nonperturbative level. 

Considerable progress has been achieved recently in applications of the operator 
product expansion to the calculation of the heavy meson decay width. In particular, 
it was shown that there are no 0(Aqcd/Mq) corrections to the Born term and first 
nonperturbative corrections 0(Aq CD /Mq) were calculated explicitly f|, 0. However, 
these contributions are all nonvalence and so have nothing to do with lifetime differ- 
ences. They are important however for the calculation of the absolute decay rates. 
There appeared a number of papers where these results (neglecting the four-fermion 
operator contributions) were applied to determine the values of the quark masses, M c 
and M b , and \V c b\ ||, @, [fjlj. However, because the four-fermion operator contribu- 
tions are of crucial importance for the D mesons and lead to t(D + )/t(D°) ~ 2.5, the 
real accuracy of the above results remained unclear. 1 



3. The heavy quark masses 

As most of calculations with heavy quarks are highly sensitive to precise values 
of their masses, it is of prime importance to know these as precisely as possible. 

Below the pole masses, M c and Mb, are used for the charm and bottom quarks 
as most convenient gauge and renormalization scheme independent quantities. Be- 
sides, the most convenient expansion parameter of the Heavy Quark Effective Theory 
(HQET) ||12|| is just the quark pole mass. On the other hand, because the pole mass 
receives contributions from the infrared region, the perturbative series connecting it 



with the current mass (say, MS'-mass) is divergent due to renormalon effects [13 

2 



+ 



Because the series in Eq.(l) is divergent, the result contains an ambiguity of order 
O(Aqcd) and requires for a concrete definition. It is possible, for instance, to cut 
out the series at the optimal number of terms: n ~ 2n / (b a s ). In what follows we 
choose the definition of the quark pole mass which looks most natural, i.e. the real 
part of the analytically continued Borel transform (Re[BT] ) of the above series (see 

1 In fact, the authors of ^ renounced their results on M c and Mf, and used M c = 1.35 — 
1.40 GeV, M h = 4.8 Gey in their later articles O. 
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i.e. PI). This was calculated in |]15| (in the improved leading order approximation, 
ILO, i.e. by replacing constant a s by running a s in one loop corrections), so that 
the result for Kq 1 ^ 1 in Eq.(l) is known. 

To illustrate the characteristic numbers, let us take (see below): M c = (1.65 ± 
0.05) GeV, M b = (5.04 ± 0.05) GeV. Then (a s is in the MS-scheme: a s (M c ) = 
0.310, a s (M b ) = 0.204) : - 1 ^ - (m) 



1.75, k£"> =2.10 |TJ: 
4« S (M C 



1 + 



1 + 1.04 + 



7T 



4q s (M c ) 
1+ — ^^1.75 

3 7T 



1.23 M r 



(2) 



M 6 = M fe (l + ^^(l + 0.62 + 
y 3 7r 



M h 



1 + — ii_^ 2 .10 

3 7T 



1.18 



(3) 



* eu + X) 
1.27 GeV 
5.04 GeV 



Therefore, the value M c = 1.65 GeV found below from the data on T(D 
leads to: M c = 1.34 GeV, which is ~ 70 MeV above the value M c = 
obtained long ago [R]] from the charmonium sum rules. The value M b - 
obtained below from M c = 1.65 GeV and HQET mass formulae leads to: M b = 
4.27 GeV, in good agreement with the value M& = 4.25 GeV obtained in jni from 
the upsilonium sum rules. 

It seems natural to expect that the accuracy of M b obtained from the upsilonium 
sum rules is better in comparison with M c obtained from the charmonium. So, 
supposing the value M b = 4.27 GeV is sufficiently accurate (say, within a several 
tens MeV), we can then turn around the chain of reasoning. Starting from this value 
we obtain M b = 5.04 GeV from Eq.(3), 2 and M c = 1.65 GeV from the HQET mass 
formulae, with each step having a several tens MeV typical accuracy. 



4. D + — > Iv + X. Determination of M c . 

On account of radiative, 0(Aq CD /M^) and four-fermion operator (0(Aq CD /M^) 
corrections, the expression for T s i can be represented in the form ||: 

orn Irad \ %o N c 

2 It seems, this value of M b is stable even on account of 0(a 2 s ) corrections. On the one hand, there 
appeared (preliminary) results ]l7t on the 0(a 2 )-corrections to the sum rules . They are positive and 
enter with significant coefficients. So, they will increase slightly the above value = 4.25 GeV. 
On the other hand, the genuine 0(a 2 )-correction in Eq.(3) is known and is negative, so that the 
transition coefficient between M b and Mf, is slightly less than 1.18. These two effects tend to 
compensate each other, so that the value = 5.04 GeV has real chanses to be sufficiently precise. 



1 - 



2 ZX rl G 

N c z M? 



+ 41 ± A 



(4) 
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Born 



G\Ml 

192 7T 3 



< D\ cc\D > 



2M, 



D 



/4-(p 2 ) ' 

2M? 



(5) 



(P 2 ) 



< D\c(iT» 2 c\D > 



2M 



D 



/4 



< D\c^g s a av GavC\D > 



2M, 



D 



(6) 



l rad 



2 a s (M c 



71 



7o«. 



(t») 
c 



(7) 



where z Q , Z\, f Q are known functions of m s /M c , 4 w) is the analog of /4 m > in Eq.(l) 
and 8\f pt in Eq.(4) is the contribution of the four-fermion operators. 
We use below 3 : 4 

m s ~ 150 MeV; (p 2 ) ~ 0.25 GeV 2 ; /i 2 G ~ \m c (M d * - M D ) ~ 0.35 GW 2 . (9) 

The value of M c can be determined from a comparison of Eq.(4) with the experi- 
mental value: T sl (D + -> Zz/ + X) = (1.08 ±0.06) • 10" 13 Gel/ . Because the dependence 
of T s i on M c is highly nonlinear, it is more convenient to proceed in an opposite way. 
Namely, let us show that Eq.(4) reproduces the experimental value at M c ~ 1.65 GeV. 
We have: 

a s (M 2 ) ~ 0.310, / ~3.25, JV C ~ 1.02 , ^tom ~ 2.80 • 10~ 13 GeV , (10) 



r^ + ) * rl rn j rod {o.7i + 41 ± a} • (ii) 

As for k^"\ the two loop radiative correction 0(b o a 2 ) was calculated in [21]: 



L rad 



2a s (M c 



71 



■fo (1 + 1.03 + 



(12) 



3 Our definition of (p 2 ) (and, analogously, A and others) differs from the cut off dependent (p 2 ) A 
used by I.Bigi et.al. If one represents this last (at /i ^> Kqcd) as: 



<p 2 ),=G,% 2 

7T 



CiA 2 



QCD 



C 2 



A 3 - 



QCD 



(8) 



then our (p 2 ) corresponds to the cut off independent term CiAq CD . We prefer this definition because 
it selects the universal number, while all /i-dependent terms will be canceled finally in observable 
quantities like decay widths. On the other hand, with our definition some useful inequalities like 
(p 2 )^ > (a*g)p wm l> e l° st > m general. 

4 From our viewpoint, the value: (p 2 ) ~ 0.5 — 0.6 GeV 2 obtained in |Q is overestimated. Let us 
recall jl9| that the mean value of the vacuum quark 4-momentum squared is: (— k 2 ) Q = 4/3 (k 2 ) Q ~ 
OA GeV 2 , and the quarks in the pion have their momenta somewhat less on the average than the 
vacuum quarks pcfl . Let us point out also that (p 2 ) enters here to N c only and plays no essential 
role. 
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Comparing with Eq.(2), it is seen that both series follow the same pattern. Besides, 
because the leading renormalon is the same in both series, it seems clear that k^"' 
will be close to k^. 5 Therefore, we estimate: k^ ' ~ 1.75, so that I ra d — 0.626. 

The main contributions to 8\f pt in Eq.(4) originate from the figs.la,b diagrams 
which show the matrix elements of the four-fermion operators in the effective La- 
grangian |TJ: 

{ 5 ^} fi9S . la , lb ^ {- 8 - 5% + i^m fl9 , b - -12% . (is) 

There are also other small contributions to S^L, the typical one is from the four- 
fermion operator hidden in the Born operator cp pc, fig.lc 6 : 

i S '>*}f t ,u = < WW 1 + ^T c(0) ' J "^ D >= 1% ' < 14 > 

Estimates show that the typical value of next corrections (denoted by ±A in Eq.(4)) 
is a few per cent. To be safe, we take: A = 6%. So : 

T s i(D + ) ~ 0.375 Y% orn (1 ± 10%) ~ (1.05 ± 0.10) ■ 10~ 13 GeV , (15) 

in agreement with data. Because the decay width is highly sensitive to the precise 
value of M c , 7 this later is tightly constrained: 8 

M c — (1.65 ± 0.05) GeV . (16) 

We were going in this section into calculation details as the precise value of M c is 
the base of all further calculations, and because the above value of M c is essentially 
higher than the commonly accepted at present values M c = 1.35 — 1.45 GeV. 



5. Mass formulae: Mb and A. 

The HQET mass formulae look as ( M = (M P + 3M V ) /4 ; 



M b - M c = M B - Mn + 



(p!>_ (e!> 

■>.\I, 2Mb 



+ 



'A 3 r 



QCD 

M? 



(17) 



M B = M b + A + (P 1 } 7/ G +o{ , (18) 



2Mb ' V M b 



5 This is supported also by the examples considered in |jl5| . 

6 Our result here differs by the factor 1/2 from those obtained in (22). 

7 For instance, with M c — 1.55 GeV the calculated T s i(D) will be more than 40% smaller. 

8 Really, higher order corrections are arranged in Eq.(4) in such a way to obtain a minimal possible 
value of M c . Therefore, M c — 1.65 GeV is rather a lower bound. 
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The expected accuracy of Eq.(17) is a several tens MeV, 9 and even better in Eq.(18). 
We obtain therefore (with the accuracy ±50 MeV, determined mainly by the uncer- 
tainty in M c : 

M b = 5.04 GeV, A = 250 MeV . (19) 

These results differ significantly from the widely used values: M b = 4.8 GeV, A = 
500 MeV. 

6. and fe . 

The knowledge of precise values of the decay constants fo and fs is of crucial 
importance for many calculations with the D and B mesons (analogously to f n for 
the pion). The calculated values of fo, fs are highly sensitive to the precise values 
of M c , M b (more precisely, to M D — M c ~ M B — M b ~ A), and increase quickly with 
increasing A. 

The QCD sum rules for the chiral current correlator which is a difference of the 
pseudoscalar and scalar current correlators, posses the advantage of being protected, 
in the chiral limit, against pure perturbative contributions (which are poorly con- 
trolable in separate correlators), and having no significant loop corrections to non- 
perturbative contributions. Using them and the above values of M c , M b , we obtain 

& 

/d(M c )~ 165 MeV, f B (M b )~120MeV, (20) 

(with the expected accuracy about 10%, which is always difficult to estimate reliably 
when dealing with the QCD sum rules). 

While the above value of fo is only slightly below the widely accepted at present 
value ~ 180 MeV, the above value of Jb is much smaller than the widely used values 
~ 180 - 200 MeV. 10 

7. Difficulties with the naive lifetime estimates. 

On account of the Born term, the leading radiative and 0(Aq CD /M 2 ) corrections 
(all nonvalence), the D meson hadronic width is {V Bom — G 2 F M 2 /6An 3 ): 11 

Kl - [ 1-50 ) rad Z N c [ 1.07 } fla ~ 1.54 Y B orn ■ (21) 

9 Supposing that (p 2 ) is within the reasonable interval, say: 0.25 ± 0.05 GeV 2 . 
10 As for the QCD sum rule calculations, the large value of originates mainly from using small 
values of Mi, (i.e. A ~ 500 — 600AfeK). As for the lattice calculations, the predictions for fs 
decrease with time, starting from ~ 250 — 300 MeV 1-2 years ago and reaching now 148 ± 20 MeV 
in the latest paper f23|j . 

11 Unlike the semileptonic width, there are two ^Lq/M 2 corrections, each one ~ 30% but they cancel 
strongly each other ||] . 
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At the level 0(A 3 QCD /M 3 ) there appear first valence (and additional nonvalence) 
contributions to the decay widths originating from the four-fermion operators. The 
effective Lagrangian (normalized at /i 2 = 0.5 GeV 2 ) has the form M: 12 

Kss-^ {^g»»L% + \ 2 t, u (\)L; u + L PNV + ---} , (22) 
L% = (-1.1 ( cr M d) ( dT v c ) + 4.0 ( cT^ d ) ( dr„¥- c ) 1 , (23) 



L; v = S.3(cT^u)(uT u ^c), t, v {X) = U^-g au ) , (21) 



where A (or A ) is the total 4-momentum of the integrated quark pair. It can be read 
off from each diagram in fig. 2 and differ from P c by the spectator quark momenta. 
The term Lp^y is nonvalence and originates from the diagram in fig. la. 

Let us try now to obtain the estimate of < D + \L e ff\D + > using the factorization 
approximation and neglecting the spectator quark momenta in comparison with M c 



24|. Then, L u and the second term in L d give zero contributions and 



ATf actor (D 



1.1 • 16tt 2 



MP 



Born 



1.50 V Born . (25) 



The term L PNV in Eq.(22) adds ~ — 0.15 T Born [0, so that we obtain for the hadronic 
width of D+: 

F nl (D + ) ~ (1.54 - 1.50 - 0.15) T Born ~ -0.11 T Born , (26) 

which does not make much sense. It is clear that the above approximations are too 
poor and some estimates are essentially wrong. 

8. Non-factorizable contributions, A and A. 

It is seen from Eqs.(22)-(24) that (internally) coloured operators whose matrix 
elements are zero in the factorization approximation, enter L e ff with much larger 
coefficients. So, even if their matrix element are suppressed, they may be of impor- 
tance. This is really the case. For instance, the most important weak annihilation 
non-factorizable contributions shown in fig. 2 were estimated using the QCD sum rules 



<D{p)\{cT^q){qT v ^c)\D(p)> iiM , (27) 



12 



Only the most important terms are shown. 
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<B(p)\{bT^q){qT^b)\B(p)> }ig ^ , (28) 

(the contribution due to the fig. 3 diagram is smaller, the fig. 2c contribution is negli- 
gibly small). Comparing with the factorizable contribution: 

<D(p)\(cV,q)(qV v c)\D{p)>^flMl(^^ , (29) 

we see that the factorization approximation works very well, even for the D mesons. 

But there is one more factor which suppresses heavily the factorizable contribu- 
tions in the D meson case: the phase space of the integrated quark pair is much smaller 
for the crossed contributions (fig. 2b) in comparison with the direct ones (fig. 2a). I.e., 
the characteristic value of A 2 in Eq.(22) is much larger than those of A . Really : 
A ~ P c + k = P D , while A ~ Pd — k\ — k 2 , where ki, k 2 are the spectator quark 
momenta. It can be shown that the momentum fraction carried by the spectator 

quark is: (x) ~ 2 (l — Mq/Mp^j /3 (Mp is the meson mass), and is ~ 15% for 

2 

the D meson and ~ 6% for the B one. So, the simplest estimate of (A ) looks as: 
(A ) /Mp ~ [1 — 2(x)] 2 and is ~ 0.5 for the D meson and ~ 0.8 for the B one. More 
accurate calculations give even stronger effect for the D meson: 

(A 2 ) D ~ 0.35 Ml, (X 2 ) B ~ 0.80 M 2 B . (30) 

As a result of all the above effects, the non-factorizable contributions become com- 
parable with the factorizable one. 

9. £) ±,0 ' s decay widths. 

Accounting for all the above described contributions, one obtains (in units 
10~ 13 GeV , the experimental values are given in brackets, Tbotu = 8.4): 

r tat (D + ) ~ T B orn [1-54 - 0.27 - 0.62 - 0.21 + (0.26)^] = 

= 0.70r flom = 5.9 {6.2}, (31) 

T tot {D°) ~T Bo m{ 1-54 - 0.27 + 0.33 + (0.26) iept ] = 

= 1.86 V Bo rn = 15.6 {15.8}, (32) 

T tot {D s ) ~T Bo m [1-54 - 0.27 - 0.07+ (0.29) iept ] = 

= 1.49 T Born = 12.5 {13.8}. (33) 

In Eqs.(31)-(33): the origin of "1.54" is explained in Eq.(21), and the term "-0.27" 
is the nonvalence contribution from figs. la, lb diagrams ("-0.15" and "-0.12" respec- 
tively). The valence terms: "-0.62" in Eq.(31) is the factorizable interference contri- 
bution (fig. 2b without gluons); "-0.21" in Eq.(31) is the summary non-factorizable 
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contribution of fig. 2b and fig. 3 diagrams; "0.33" in Eq.(32) is the non-factorizable 
direct annihilation contribution, fig.2a; the small term "-0.07" in Eq.(33) originates 
from the Cabibbo-suppressed contributions and the non-factorizable annihilation con- 
tributions. 

The expected accuracy of T tot (Di) calculations is not high: ~ 20%. 13 It is 
sufficient however to keep all the main contributions well under control, and is much 
better in comparison with all previous estimates which were qualitative rather than 
quantitative (see sect. 7). 



10. B ±,0 ' s lifetimes and B — B mixing. 

The u(d + s)— part of the B d width can be represented in the form: 

(2C 2 + + C 2 _] 



J- (ud) (B d ) ~ T B 

orn 



ud j(ud) q 
^o 1 rad ^''power i 



(34) 



where: (2C+ + C 2 )/3 ~ 1.127 is the renormalization factor due to the evolution 
from M w to M b [p5| , Q power ~ 0.99 is the summary effect of all power corrections 

rod — 1-03 p6[ describes all other 
0.13, l£J ~ 1.3 |1, ||. On the 



^ud 



0.46 is the phase space factor, and 1 



0(a s ) corrections. For the c(s + d) part: z c Q & 
whole (see sect. 11 for the semileptonic decays and V cb , T Born = GpM^\V c b\ 2 /64it 3 



3.93 ■ 10- u GeV , r Bc 

Ttot(Bd) — r Born 

~ 1.00 -T Bom - 



1.67 ps, the data see in f28 



0.53 (nd) + 0.19 



(cs) 



3.93 • 10~ 13 GeV 



2 • 0.1U (eu+uu) + 0.028 M + 0.02 (b ^ u) j 

(B d ) = (4.10 ±0.18) • 10" 13 GeV ,(35) 



-pexper 
1 tot 



Br {eu) ~ 11.4% , Br {Tv) ~ 2.8% , Br {cs) ~ 19% . (36) 



All these results agree with data 14 p8| , although Br^ is slightly above the exper- 
imental value: Brf^ er) = 0.11 ± 0.06. 15 

As for the lifetime differences, the qualitative picture remains the same as for the 
D mesons but, of course, all the effects are much smaller. Let us denote by r o the 
common width of all B ±,0,s mesons which they have on neglect of the four-fermion 
operator contributions (and SU(3) breaking). Then the valence contributions look as: 

) ~ ( -2.1% ) + ( -1.3% ) ^ -3.4% , (37) 

J- o 



13 Besides, the SU(3) symmetry breaking corrections are not accounted for in calculations with D s 

14 One can expect also that higher loop corrections will increase slightly the hadronic width. 

15 Let us emphasize that Bri cs \ will be essentially larger (> 0.30) for the quark masses M c i 



'(cs) 

1.4 GeV, Mb ~ 4.8 GeV, and this will constitute a real difficulty. 
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where (—2.1%) is the factorizable interference contribution (fig. 2b without gluons), 
and (—1.3%) is the summary non-factorizable contribution of fig. 2b and fig. 3 dia- 
grams. 

The Bd and B s mesons (neglecting SU(3) breaking) receive only the non-factorizable 
contributions from fig. 2a and fig. 3 diagrams: 

(38) 

J- o J- o 

Finally, the four-fermion operators give also the non-valence contributions, both 
factorizable: ~ (—0.6%) from fig. la diagram, and non- factorizable: ~ (—0.1%) from 
fig. lb one. 

On the whole, the lifetime difference of the Bd and B ± mesons is [|TJ : 

T ™-£ B±) - 4%, (39) 

while Ar of Bd and B s is zero within the available accuracy. 

The above described non-factorizable contributions, figs. 2, 3, determine also the 
corrections to the factorization approximation for the B — £>-mixing, and appear to 
be surprisingly large here [|IJ: ~ —18%. As a result, one obtains for the S s -mesons 
(neglecting SU(3)-breaking): 

(40) 

and for the "bag factor" of the Bd — Bd mass mixing: 

< B d \(bT u d)(bT u d)\B d >,= Mb = y 2 B {M b )M 2 B B B {M b ) , 

B B (M b ) ~ (1 - 18%) = 0.82 . (41) 



11. B —>■ ev + X c : determination of \V C 



cb\ 



The calculation proceeds in analogy with those in sect. 4. The total effect of 
0(AQ CD /M b ) power corrections is much smaller now (~ —4%), and 5^ pt is negligibly 



small (~ -2 ■ 10" 
(f ~ 2.46 here): 



Most important are radiative corrections which look as |f2l} 



'rad 



1 _?^) /o _ 1 . 68( , o f££« 

3 7T V 7T 



2a s (M b ) , . 
1_ _!L_^/ o 1 + 0.6 + 

O 7T 



2 a s (M b ) {w) 
3 Z Jo^b 



71 



(42) 
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Comparing with Eq.(3) it is seen that (analogously to the charm case) two series follow 
the same pattern. So, we estimate: /4 ~ 2.1, I ra d — 0.776 . Therefore (zf ~ 0.46): 



r(B -> e i/ + X c ) ~ 0.44 r^ orn J rad ~ 4.5 • 10" 14 GeV 



v cb 1 



0.042 



(43) 



T si = gllMM ^ i3.i . io- 14 Gey 
Sorn 192 vr 3 



c6 



0.042 



(44) 



As for the data, we take [2S1: 



r B = ( 1.60 ± 0.07 ) ps , Br[B -> ez/ + X) = (11.0 ± 0.5)% 



(45) 



and obtain: 



Va \ = (42 ±1 ) • 10" 



10" 3 


"5r(5 -> ez/ + X)" 


1/2 


T.6ps" 


1/2 


11.0% 




- Tg - 





(46) 



The error bars in Eq.(46) were estimated by varying: 1.60 < M c < 1.70 GeV, 
3.37 < (M b - M c ) < 3.41 GeV . 16 

For the B -> rz/ + X decays: / ~ 2.0, ^ ~ 0.105, J rad ~ 0.83, and so: 



Br ( - ] ~ 0.25. 



(47) 



12. The unitarity triangle. 

1) . Using: V cb = A A 2 , A ~ 0.221 and (see Eq.(46)) V cb ~ 4.2 ■ 10~ 2 , one has: 

A - 0.86. (48) 



2) . The B° — B° mass difference is given by the well known formula (see i.e. [p9|), 



and using ( see Eqs.(20), (41) ): f B (M b ) = 120 Me V B B (M b ) = 0.82, and |28 



Mf°' e = 175 GeU, r(B d ) = 1.6 ps, (49) 

one obtains : 

x d = — = (0.78 ± 0.06) ~ 4.5 ■ 10 3 \V td \ 2 , (50) 



16 For comparison, using M c = 1.4 GeV, Mfc = 4.8 GeV and proceeding in the same way one obtains 
44.5 instead of 42.0 in Eq.(46), and it seems it will be difficult to reconcile this value with the data 
on T(B -> D*ev). 
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Wt 



td\ 



P? + v 2 



1/2 



1.3-10" 



(51) 



2.0. 



(52) 



3) . Using the above given parameters and Bk — 0.82 from the lattice calculations 
(J, the CP- violating part of the K° — K° mixing can be written in the form (see 
i.e. P|): 



e ~W4 e ^ . 10 3 _ g fe K v ( L36 - p ) = 2.26 , 



7/(1.36-/o) ~ 0.35. 
Therefore, we obtain from Eqs.(52), (53): 

p~-0.40, T] ~ 0.20, 5 = arctg(-) ~ 0.85 tt . 

P 



V, 



ub 



cb 



0.10 



(53) 



(54) 



sin 2a ~ 0.60, sin 2(3 ~ 0.28, sin 27 ~ -0.80 . 



(55) 



The unitarity triangle is shown in fig. 4 . With the above parameters the CP-violating 
asymmetry in the B° d — > ^ decay is: 



|^(S d °^*K s ) 



^-^sin 2/3 ~ 0.14. 
1 + x^ 



(56) 



13. Summary 

It is seen from all the above presented calculations that a self-consistent picture 
emerges which agrees with a large number of various experimental data, and allows to 
obtain a number of important predictions which can be checked in future experiments. 

One of our main concerns was to calculate reliably the four-fermion operator 
contributions. These are of crucial importance for explaining the pattern of the 
D±,o,s lifetimes, the B ±,0,s lifetime differences and the B — B mixing. There is a 
clear reason explaining the importance of these four-fermion operator contributions, 
although they are formally only 0(Aq CD /Mq) corrections: they are the first who 
gain the large numerical factor ~ 167r 2 (see Eq.(25)) due to the two-particle phase 
space, in comparison with the three-particle one for the Born term and (non- valence) 
0(Aq CD /Mq) corrections. It is clear that this enhancement factor operates one time 
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only, so that all other 0(Aq CD /Mq) and higher order corrections are naturally small 
(see Eq.(14)) and have no much chances to be of real importance. 17 

The calculated values of the c and b quark pole masses, M c ~ 1.65 GeV, Mj, ~ 
5.04 GeV, appeared to be significantly larger the widely accepted at present values: 
M c ~ 1.35 — 1.45 GeV, Mb ~ 4.8 GeV. This difference is of great importance, as 
most of calculations with heavy quarks are highly sensitive to precise values of their 
masses. In particular, the calculated decay constant fs — 120 MeV appeared to be 
much smaller the widely accepted at present value fs — 180 — 200 MeV. This dif- 
ference leads, in its turn, to essentially different predictions for the B° — B ± lifetime 
difference, B — B° mixing and the unitarity triangle parameters. Just because the 
above obtained values of M c , Mb, A, etc, look highly non-standard at present, 
we described above the calculations of many experimentally measured quantities to 
show there is no disagreement with data. Besides, a number of concrete predictions 
is described which can be checked in future experiments (see Tables). 
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17 The applicability of the standard operator expansions to the calculation of T(D — > ev + X) has 
been questioned in (22) on the only ground that the authors don't believe the c-quark pole mass 
may be as large as M c ~ 1.65 GeV. They insist it can not exceed ~ 1.4 GeV. Let us emphasize that 
the calculated value of T(D — > ev + X) (see Eq.(15)) will decrease ~ 3 times for M c = lAGeV. So, 
there should be huge additional contributions which dominate the semileptonic width and remain 
invisible within the standard operator expansion. No one reliable argument is presented however in 
p3] to justify M c ~ 1.4 GeV, and no one missed contribution is shown which is of great importance. 
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Table 1: INPUT 



T(D -> ev + X) 


(1.08 ±0.06) • 10~ i3 GeV 


r{B d ) 


(1.60 ±0.07) ps 


Br(B -> ez/ + X) 


(11.0 ±0.5)% 




0.78 ±0.06 




2.26 • 10~ 3 




0.82 ±0.05 


Mf°' e 


175 GeV 


A = |Vy 


0.221 




0.118 


\f B B B \ B J\f 2 B B B \ Bd 


1.3 



Table 2: OUTPUT 



M pole 


1.65 GeV 


M pole 


5.04 GeV 


fn(M c ) 


165 MeV 


f B (M b ) 


120 MeV 


B B (M b ) 


0.82 




T(B d ) - T(B±) 




4% 






AW 


8% 




T{B s d hort ) - T(B l ° ng ) 


m) 


0.6% 


Am s /Am d 


14 


\Vcb\ 


4.2 • 10~ 2 


\Vtd\ 


1.3 • 10~ 2 


\v ub /v cb \ 


0.10 


\A(B^^K S ) 




0.14 


\e'/e\ • 10 4 


a few units 



Figure captions 

Fig.l Non- valence factorizable (a,c) and non-factorizable (b) contributions to 
matrix elements. 

Fig. 2 Valence direct (a) and cross (b) annihilation non-factorizable contribu- 
tions to matrix elements of coloured operators; (c) the same for colourless operators. 
Fig.3 Valence non-factorizable contribution (plus the mirror diagram). 
Fig.4 The unitarity triangle: p ~ —0.40 , r] ~ 0.20 . 
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